Interface mediated interactions between particles — a geometrical approach 
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Particles bound to an interface interact because they deform its shape. The stresses that result 
are fully encoded in the geometry and described by a divergence-free surface stress tensor. This 
stress tensor can be used to express the force on a particle as a line integral along any conveniently 
chosen closed contour that surrounds the particle. The resulting expression is exact (i. e., free of 
any "smallness" assumptions) and independent of the chosen surface parametrization. Additional 
surface degrees of freedom, such as vector fields describing lipid tilt, are readily included in this 
formalism. As an illustration, we derive the exact force for several important surface Hamiltonians 
in various symmetric two-particle configurations in terms of the midplane geometry; its sign is 
evident in certain interesting limits. Specializing to the linear regime, where the shape can be 
analytically determined, these general expressions yield force-distance relations, several of which 
have originally been derived by using an energy based approach. 

PACS numbers: 87.16.Dg, 68.03. Cd, 02.40.Hw 
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I. INTRODUCTION 

The interaction between spatially separated objects is 
mediated by the disturbance of the region that surrounds 
them, described by a field. In electromagnetic theory for 
example the interaction between charged particles is de- 
scribed by the Maxwell field equations. Since they are lin- 
ear, interactions add. However, more often than not, the 
field equations are nonlinear as for example in the case of 
General Relativity: even though the energy-momentum 
tensor couples linearly to the curvature, the latter de- 
pends in a nonlinear way on the spacetime metric and 
its derivatives [ij, |2|, |3( • The source of the nonlinearity 
lies in the geometric nature of the problem. Not only do 
interactions fail to add up, even the humble two particle 
problem poses challenges. 

"Effective" interactions between macroscopic degrees 
of freedom arise in statistical physics when a partial 
trace is performed in the partition function over unob- 
served microscopic degrees of freedom |j, |5( . The Boltz- 
mann factor invariably renders these interactions non- 
linear. This time, the source of the nonlinearity is the 
entropy hidden in the degrees of freedom that have been 
traced out. For example, the effective interaction be- 
tween charged colloids in salty water is described (at a 
mean-field level) by the Poisson-Boltzmann equation Q . 

In this paper we will discuss a classical example which 
belongs to the class of effective interactions, while ow- 
ing its nonlinearity to its geometric origin: the inter- 
action between particles mediated by the deformation 
of a surface to which they are bound. This prob- 
lem includes the capillary interactions between particles 
bound to liquid- fluid interfaces |7], la, |9|, or the mem- 
brane mediated interactions between colloids or proteins 
adhering to or embedded in lipid bilayer membranes 



HlliimEEQHm- To approach the problem 
we require two pieces of information. First, how does the 
energy of the surface depend on its shape, or in other 
words, what is the "surface Hamiltonian" ? Second, how 
does a bound particle locally deform the surface? With 
this information, we may (in principle) deduce the equi- 
librium shape minimizing the energy of the surface for 
any given placement of the bound particles. Knowing 
the shape, the energy can be determined by integration, 
and the forces it implies follow by differentiating with 
respect to appropriate placement variables. In general, 
however, the ground state of the surface is a solution of a 
nonlinear field equation ( "the shape equation" ) , thereby 
thwarting progress by this route at a very early stage. 

Sometimes the linearization of a nonlinear theory is 
adequate. Just as one recovers Newtonian gravitation 
as the weak- field limit of General Relativity 0, S > ° r 
Debye-Hiickel theory as the weak-field limit of Poisson- 
Boltzmann theory 6j, a linear theory for surface medi- 
ated interactions is useful for certain simple geometries, 
notably weakly perturbed flat surfaces. At this level, 
the approach to interactions based on energy becomes 
tractable. Yet, linearization is also often inadequate. The 
full theory may display qualitatively new effects which 
are absent in the linearized theory: strong gravitational 
fields give us black holes |J H @; the bare charge of 
a highly charged colloid gets strongly renormalized by 
counterion condensation [T7j ]. 

There is, however, an alternative approach to interac- 
tions which was outlined in 18] . By relating the inter- 
action between particles to the equilibrium geometry of 
the surface, a host of exact nonlinear results is provided. 
The link is formed by the surface stress tensor, and it can 
be established without solving the shape equation. Once 
we know the Hamiltonian, we can express the stress at 
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any point in terms of the local geometry - covariantly 
and without any approximation. We will briefly revisit 
the essentials of this construction in Sec. [H] Knowing 
the stresses, the force on a particle is then determined 
by a line integral of the stress tensor along any surface 
contour enclosing the particle, as we will show in Sec. II VI 

Such results might, at first sight, appear somewhat for- 
mal: without the equilibrium surface shape, they cannot 
be translated into hard numbers. However, the close link 
between the force and the geometry, combined with a 
very general knowledge one has about the surface shape 
(e.g., its symmetry) will turn out to provide valuable 
qualitative insight into the nature of the interaction (e. g., 
its sign). Even on a completely practical level, this ap- 
proach scores points against the traditional approach in- 
volving energy, providing a significantly more efficient 
way to extract forces from the surface shape determined 
numerically (in whatever way). 

We will illustrate this approach with a selection of ex- 
amples involving different symmetries and surface Hamil- 
tonians. In Sec. lIIII we demonstrate how its scope extends 
in a very natural way to include internal degrees of free- 
dom on the membrane - in particular: a vector order 
parameter which has for instance been used to describe 

lipid tilt m m m, m s 113. To make contact with 

the energy based approach in the literature, and also in 
order to link the formalism to a more familiar setting, 
we specialize in Sec. to a Monge parametrization and 
its linearization. This will permit us in Sec. IVII to de- 
rive force-distance curves for interactions mediated by 
surface tension, membrane curvature, and lipid tilt. Var- 
ious well-known linear results 0, H, UM then follow very 
naturally using the stress tensor approach. 



II. ENERGY FROM GEOMETRY 



In this paper we want to study the physics of inter- 
faces, which are characterized by a reparamctrization in- 
variant surface Hamiltonian. The appropriate language 
for this is differential geometry, and in this section we 
will outline how physical questions can be formulated 
very efficiently in this language. We first summarize the 
necessary mathematical basics and introduce our nota- 
tion (the reader will find more background material in 
Refs. |23). We then define the class of Hamiltonians we 
will be considering. The corresponding Euler-Lagrange 
equations for the interface degrees of freedom will be cast 
as a conservation law. The most direct way to do this is 
to implement all geometrical constraints using Lagrange 
mulipliers; not only does this approach provide a quick 
derivation of the shape equation, it also provides a trans- 
parent physical identification of the surface stresses. 



A. Differential Geometry and Notation 

We consider a two-dimensional surface £ embedded in 
three-dimensional Euclidean space R 3 , which is described 
locally by its position X^ 1 ,^ 2 ) e R 3 , where the £ a are 
a suitable set of local coordinates on the surface. The 
embedding functions X induce two geometrical tensors 
which completely describe the surface: the metric g a b 
and the extrinsic curvature K ab , defined by 

g a b = e a ■ e b and (la) 
K a b = e a ■ d b n , (lb) 

where a, b E {1,2}. The local coordinate frame formed 
by the tangent vectors e.\ and extended by the normal 
vector n forms a local basis of M 3 : 

e a = dX/dC = d a X , (2a) 
e a n = , (2b) 
n 2 = 1 . (2c) 

Note that unlike n, the e a are generally not normalized. 

In the following, V„ denotes the metric-compatible co- 
variant derivative j2g and A = V a V a the correspond- 
ing Laplacian. Surface indices are raised with the in- 
verse metric g ab . The trace of the extrinsic curvature, 
K = 9 ab Kab, is twice the mean curvature. Using the 
above sign conventions, a sphere of radius a with out- 
ward pointing unit normal has a positive K = 2 /a. 

The intrinsic and extrinsic geometries are related by 
the Gauss-Codazzi-Mainardi equations 

V a K bc - V b K ac = , (3a) 
K ac Kbd — K ad K bc — R abc d , (3b) 

where R a bcd is the Riemann tensor constructed with the 
metric; its contraction over the first and third index is 
the Ricci tensor, Rbd = g ac R a bcd, whose further con- 
traction gives the Ricci scalar curvature, R = g bd Rbd- 
From Eqn. (|3b|) we see that the latter satisfies R = 
K 2 — K ab K ab . In particular in two dimensions we have 
that R — 2 Kq, where Kq = det(i^) is the Gaussian 
curvature (Gauss' Theorema Egregium (25|). 

B. Surface energy and its variation 

We consider surfaces such as lipid membranes and soap 
films, characterized by the property that the associated 
energy is completely determined by the surface geometry 
and described by a Hamiltonian which is an integral of a 
local Hamiltonian density Ti over the surface: 

H[X] = [ dA Ti(g ab , K ab , V a if fcc , . . .) . (4) 

The infinitesimal area element is dA — y^d 2 ^, where 
g = det(g ab ) is the determinant of the metric. The den- 
sity Ti depends only on scalars constructed from local 
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surface tensors: the metric, the extrinsic curvature, and 
its covariant derivatives. In order to find the equilibrium 
(i. e., energy minimizing) shape, one is interested in how 
H responds to a deformation of the surface described by 
a change in the embedding functions, X — > X + SX. 
The straightforward (but tedious) way to do this is to 
track the course of the deformation on X through g ab , 
yfg, K ab , and any appearing covariant derivatives using 
the structural relationships and 101. 

Alternatively, one can treat g ab , K ab , e a and n 
as independent variables, enforcing the structural re- 
lations Q and J2J) using Lagrange multiplier func- 
tions 27]. One thus introduces the new functional 
H c [g ab , K ab , ...,X,e a , n, A afc , A ab , f a , X±, A„] given by 

H c = H[g ab ,K abl . . .] 
Gl f 

+ J dA \ ab (g ab - e a ■ e b ) 

ITEi 

+ 

+ 

l2bl 

+ 

He} 



+ J dAA ab (K ab -e a -d b n) 
dA f a ■ (e a - d a X) 

2b) f 

+ / dAAl(e Q -n) 
dA A„(n 2 - 1) . 



(5) 



The original Hamiltonian H is now treated as a func- 
tion of the independent variables g ab , K ab and its covari- 
ant derivatives; X ab , A ab , f a , A° and A„ are Lagrange 
multipliers fixing the constraints and (J2J). The intro- 
duction of auxiliary variables greatly simplifies the vari- 
ational problem, because now we do not have to track 
explicitly how the deformation SX propagates through 
to g ab and K ab . As we will see in the following, this ap- 
proach also provides a very simple and direct derivation 
of the shape equation in which the multiplier f a , which 
pins the tangent vectors to the surface, is identified as 
the surface stress tensor. 

Note that additional physical constraints can be en- 
forced by introducing further Lagrange multipliers (such 
as a pressure P in the case that a fixed volume is enclosed 
by the surface). 



Euler-Lagrange equations for X, e Q , n, g ab and K ab : 

V a / a = , (6a) 

f a = (A ac K b + 2X ab )e b - X a ± n , (6b) 

= (V 6 A ah + Al)e a + (2A„- A ab K ab )n , (6c) 
X ab = i T ab ; ( 6d ) 

A ab = _ n ab (g e ) 



Note that the Weingarten equations d a n = K b e b have 
been used in Eqn. (|6b(l : the Gauss equations V a e6 = 
—K ab n have been used in Eqn. JJjEJl. We have also defined 



7~L ab 

rpab 



5H 

and 

5K ab 

2 S(^gH) 

Vg bgab 



(7a) 
(7b) 



The manifestly symmetric tensor T ab is the intrinsic 
stress tensor associated with the metric g ab . If TL does 
not depend on derivatives of K ab , functional derivatives 
in the definition of Ti ab and T ab reduce to ordinary ones. 

Equation l|6a|l reveals the existence of a conservation 
law for the current f a . Using the other equations (JScJ, 
(|6d|) and (|B*e|l . it is straightforward to eliminate the La- 
grange multipliers on the right hand side of Eqn. (|6bl) 
to obtain an explicit expression for f a in terms of the 
original geometrical variables. From Eqn. (|6c|l we find 
A^ = — VftA ab because e a and n are linearly indepen- 
dent; the Eqns. (|6d|) and l|He|l determine A ah and A ab . 
Thus Eqn. 16 bp can be recast as 



f a = (T ab - H ac K b )e b - {W b H ab )n 



(8) 



Once the Hamiltonian density has been specified, 
Eqn. (jHJ) determines the conserved current f a completely 
in terms of the geometry. Several representative exam- 
ples are treated in the Appendix. 

Finally, as pointed out in Refs. HSH^i the normal pro- 
jection of V a / a is the Euler-Lagrange derivative £(TL) of 
the original Hamiltonian H which vanishes for an equilib- 
rium shape |29|. Using the Gauss equations once more, 
we obtain the remarkably succinct result 



n ■ V a / a - £ (H) = -K ab T ab + (K ac K c b - V a V b )H afc 



(9) 



C. Euler-Lagrange equations and the existence of a 
conserved current 

The Hamiltonian is invariant under translations. 
As explained in detail in Ref. |28|, Noether's theorem 
then guarantees the existence of an associated conserved 
current, which we will identify as the surface stress tensor 
in Sec. Ill Dl In order to see this, let us first work out the 



D. Identification of the stress tensor 

We will now show that f a can be identified with the 
surface stress tensor. The variation of the Hamilto- 
nian has a bulk part proportional to the Euler-Lagrange 
derivative as well as boundary terms: under a change 
in the embedding functions X — > X + SX one gets 

SH C = [dA [V a / a • SX - V s (/° • SX)] . (10) 
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FIG. 1: Surface E with 3 disjoint boundary components <9E; 
and an outer limiting boundary <9E ou t- 



Additional boundary contributions stem from the varia- 
tions with respect to n, g ab and K ab , since these terms 
do or may contain further derivatives which then need 
to be removed by partial integration. However, the one 
appearing in Eqn. i|l(J|) is the only one that is relevant 
for identifying the stress tensor: As we will see below, 
for this we are exclusively interested in translations, for 
which n, g ab and K ab remain unchanged. 

Consider, in particular, a surface region £ in equilib- 
rium (see Fig.^l: its boundary <9£ consists of n disjoint 
closed components dYn and an outer limiting boundary 
<9£out ■ Each of the dY*i is also the closed boundary of a 
surface patch S^. Under a constant translation 5X = a 
of dYn the only non-zero term is 

5H C = -a l ds l*f a = a ■ F« . (11) 

Stokes' theorem has been used to convert the surface in- 
tegral into a line integral. The vector i E = 1^ e a is the 
outward pointing unit normal to the boundary on the sur- 
face S; by construction it is tangential to S. The variable 
s measures the arc-length along <9£i. The boundary inte- 
gral is thus identified as the external force -T^xt acting on 
dYii. dotted into any infinitesimal translation, it yields 
(minus) the corresponding change in energy |30| . 

The external force -F C xt on the surface patch Eo is sim- 
ply given by — F^ t due to Newton's third law 

Fext = / ds l a f a = [ dA V a / Q , (12) 

JdY, "'So 

where I = — i s and Stokes' theorem was used again. 

Recall now that in classical elasticity theory |3l| the 
divergence of the stress tensor at any point in a strained 
material equals the external force density. Or equiva- 
lently, the stress tensor contracted with the normal vec- 
tor of a local fictitious area element yields the force per 
unit area transmitted through this area element. Com- 
paring this with Eqn. ljT2l we see that f a is indeed the 
surface analog of the stress tensor: l a f a is the force per 
unit length acting on the boundary curve due to the ac- 
tion of surface stresses. 

It proves instructive to look at the tangential and nor- 
mal projection of the stress tensor by defining 

f a = f ab e b + f a n . (13) 



Using the equations of Gauss and Weingarten |32(, the 
relation V / a = £n can then be cast in the form 

V a / a = K ab f ab + £ , (14a) 
V Q / Qb = -K b J a . (14b) 

Tangential stress acts as a source of normal stress - 
and vice versa. Both conditions hold irrespective of 
whether the Euler-Lagrange derivative £ actually van- 
ishes. In fact, Eqn. I|14al) shows that the shape equation 
£ = is equivalent to V Q / a = K ab f ab , while Eqn. (fl4E|) 
merely provides consistency conditions on the stress com- 
ponents. For instance, the Helfrich Hamiltonian H cx K 2 
yields f a cx W a K, while f ab is a quadratic in the extrin- 
sic curvature tensor (see Eqn. ljA4p ). Hence, Eqn. 114a|) 
immediately reproduces the characteristic form of the 
Euler-Lagrange derivative: AK plus a cubic in the ex- 
trinsic curvature. 



III. INTERNAL DEGREES OF FREEDOM 

So far we have restricted the discussion to Hamiltoni- 
ans which are exclusively constructed from the geometry 
of the underlying surface. However, the surface itself may 
possess internal degrees of freedom which can couple to 
each other and, more interestingly, also to the geome- 
try. The simplest example would be a scalar field <f> on 
the membrane, which could describe a local variation in 
surface tension or lipid composition, and it is readily in- 
corporated into the present formalism |33|. 

Here we will look a little more closely at the case of 
an additional tangential surface vector field m a . Such a 
field has been introduced to describe the tilt degrees of 
freedom of the molecules within a lipid bilayer, to ac- 
commodate the fact that the average orientation of the 
lipids themselves need not coincide with the local bilayer 
normal (see for instance Refs. [HIM IH El HI) 
Many additional terms for the energy emerge in the pres- 
ence of a new field m a (for a systematic classification see 
Ref . |2l| ) . However our aim here is not to treat the most 
general case. Instead, we will focus on a simple rep- 
resentative example to illustrate how easily the present 
formalism generalizes to treat such situations. 

Let us define the properly symmetrized covariant tilt- 
strain tensors M ab and F ab according to 

M ab = I(v a m 6 + V 6 m a ) , (15a) 

F ab = \7 a m b - V b m a . (15b) 

In the spirit of a harmonic theory we construct a Hamilto- 
nian density 7i m from the following quadratic invariants: 

H m = ^XM 2 + iiM ah M ab + -vF ah F ab + V{m 2 ) , (16) 

where M = g a bM ab = V a m a is the tilt divergence. The 
first two terms coincide with the lowest order intrinsic 
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terms identified by Nelson and Powers |2lJ, provided we 
restrict to unit vectors m a [34[- These terms are multi- 
plied by new elastic constants A and fx, playing the analo- 
gous role to Lame-coefficients [35J. If m 2 J=l a third term 
(also absent in usual elasticity theory |31|) occurs, the 
quadratic scalar constructed from the antisymmetrizcd 
tilt gradient; its structure is completely analogous to the 
Lagrangian in electromagnetism 36] . Finally, if the mag- 
nitude of m a is not fixed, we may also add a potential V 
depending on the square m 2 = m a m a of the vector field 
m a . Without loss of generality we assume that V(0) = 0, 
because any nonvanishing constant is more appropriately 
absorbed into the surface tension a. If V(x) is minimal 
for x — 0, then m a = will minimize the energy, but 
depending on physical conditions V may favor nonzero 
values of \m a \. This is why below the main phase transi- 
tion temperature of lipid bilayers the lipids can acquire 
a spontaneous tilt. 

This particular choice for 7i m is purely intrinsic. 
Hence, Eqn. (JSJ) shows that the corresponding mate- 
rial stress is also purely intrinsic, therefore tan- 
gential, and given by = T^e b , where = 
—2^/g~ 1 5(^/gli. m )/5g ab is the metric material stress. 
Performing the functional variation (see Appendix) we 
find 



rpab 



A(M 2 + 2m c V c M) + v(e cd \7 c m d ) 2 g ab 

- M cd M cd g ab + 2MM ab + 2m c V c M ab 

- (V c m a )(V c m b ) + (V a m c )(V b m c ) 
V{m 2 )g ab - 2V'{m 2 )m a m b , (17) 



derivative with respect to the field m a , which is given by 
SH 

£ma = = -AV a V 6 m b - (ji + v)V b V a m b 

- (/i - i/)Am a + 2V'(m 2 )m a . (19) 



In general, the equilibrium condition £ ma = implies 
Eqn. I|18b|) . For a single vector field m a the converse also 
holds so that Eqn. l|18bl) may be used in place of the 
equilibrium condition 38] . 

In equilibrium, we have not only £ raa = 0, we also 
have V a £ ma = 0. Using the commutation relations for 
covariant derivatives |39j, it is then easy to see that the 
tilt also satisfies the following equation on the surface: 

(X + 2fi)AM + fiV a {Rm a ) -2[2V"m 2 + V'M] = 0. 

(-") 

Notice that v has dropped out of this equation, which 
follows from the fact that F ab is invariant under U(l) 
gauge transformations ^(j- For small values of tilt, we 
can expand the potential as 



Vim 2 ) = -tm 2 H — iim 4 
2 4 



(21) 



In the untilted phase we can terminate this expression 
after the first term (since then t > 0). If we now restrict 
to a flat membrane (and thus R = 0) Eqn. (|2(Jfl simplifies 
to a Helmholtz equation for the tilt divergence: 



[(A + 2(i)A - t] M = 



(22) 



showing that (in lowest order) any nonzero M is (essen- 
tially) exponentially damped with a decay length of 



where £nj, = n ■ (e a x e b ) is the antisymmetric epsilon- 
tensor |32|. Notice that the metric stress tensor is 
quadratic in the tilt-strain, not linear. Unlike the stress 
tensor in elasticity theory, this tensor is not obtained as 
the derivative of the energy with respect to the strain 
but rather with respect to the metric, which leaves it 
quadratic in the strain. The formal analogy alluded to 
earlier is therefore not complete. 

Adding the material stress to the tangential geo- 
metric stress f ab , we find with the help of Eqns. 114|) the 
equilibrium conditions 





ry rpab 
V a^m 



-K ab T, 
. 



a b 



+ £ 



(18a) 
(18b) 



A + 2[i 
t 



(23) 



If t < gets us into the tilted phase, the expansion 121|) 
has to be taken one order higher, leaving instead a non- 
linear Ginzburg-Landau equation to be solved. 

We finally remark that even though the system of 
Euler-Lagrange equations ifTHl) is quite formidable, it still 
enjoys one nice nontrivial property: The material equa- 
tion l|18b(l is purely intrinsic. This is the case because 
the material stress is tangential, which itself derives from 
the fact that the material Hamiltonian is intrinsic. If we 
were to add a coupling between tilt and extrinsic curva- 
ture, such as the chiral term e ac K b 2 m a m b , this decoupling- 
would no longer hold. 



The first of these equations shows how the material de- 
grees of freedom "add" to the geometric Euler-Lagrange 
derivative £ of the geometric Hamiltonian 7i; this is the 
modified shape equation. The second equation - which 
before provided consistency conditions on the geomet- 
rical stresses - tells us that the material stress tensor 
is conserved. The equilibrium of the material degrees 
of freedom involves the vanishing of the Euler-Lagrange 



IV. FORCES BETWEEN PARTICLES 

Particles bound to an interface can exert indirect forces 
onto each other. Since these are mediated by the inter- 
face, they must be encoded in its geometry. We have seen 
that the "coding" is done by the surface stress tensor f a . 
The problem is to decode this content. 
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In this section we will solve this problem. The strong 
link between stress and geometry can be easily turned 
into exact expressions for mediated interactions. The 
method by which we obtain these results for various dif- 
ferent Hamiltonians as well as the final formulas are one 
of the major results of this paper. 



A. The stress tensor and external forces 

Consider a single simply-connected patch Eo- The ex- 
ternal force acting on it is given by Eqn. Ijl2|l . If there 
are no external [41| forces acting on Eo, the integrals 
appearing in Eqn. (|12|> will vanish; but even when F ext 
does not vanish, the stress tensor remains divergence free 
(Eqn. (jflajl ) on any part of the surface not externally 
acted upon. As a result, the contour integral appear- 
ing in Eqn. (|12f> will be independent of the particular 
closed curve so long as it continues to enclose the source 
of stress and does not encroach on any other sources. 

Observe now that in general a multi-particle configu- 
ration can be stationary only if external forces constrain 
the particle positions. These are the forces providing the 
source of stress in Eqn. (|12|l . The force F we are ulti- 
mately interested in is the force on a particle mediated by 
the interface counteracting this external force; we there- 
fore evidently have F = —F cxt . 



for the surface stress tensor associated with this Hamil- 
tonian. To facilitate the calculation of the force it is con- 
venient to introduce an orthonormal basis of tangent vec- 
tors {t, 1} adapted to the contour 3Eo: t = t a e a points 
along the integration contour and, as introduced previ- 
ously, I — l a e a points normally outward. The elements of 
the extrinsic curvature tensor with respect to this basis 
are given by 



(27a) 
(27b) 
(27c) 



We obtain for the integrand appearing in the line integral 
in Eqn. (JUJl, 





= l a l b K ab 


K \\ 


= t a t b K ab 




= l a t b K ab 



Lf a = n(l a K 



-Kl b )K-al b 
2 > 



e b -n{V±K)n , (28) 



where we have defined the normal derivative Vj_ = / a V a . 
The first term can be simplified by exploiting the com- 
pleteness of the tangent basis, — l b l c + t b t c : 



l a K ab e b 



l a K ab {l b l c + t b t c )e c 
l a l b K ab l + l a t b K ab t 
K ± l + K ± ut . 



(29) 



Since furthermore the trace K = K± + Ku , we find 



B. Force between particles on a fluid membrane 

Let us now focus on a symmetric fluid membrane, de- 
scribed by the surface Hamiltonian 



(24) 



which, up to irrelevant boundary terms, is equivalent to 
the Hamiltonians introduced by Canham |42j and Hel- 
frich |43|. Here, k is the bending rigidity and a is the 
lateral tension imposed on the boundary. For typical 
phospholipid membranes k is of the order of a few tens 
of fceT, where k-gT is the thermal energy. Values for a 
are found to be in a broad range between up to about 
10 mN/m 44j. The Hamiltonian (J23J) covers interesting 
special cases in various limits: soap films on setting k = 
and tensionless membranes on setting a = 0. Note that 
the two elastic constants provide a characteristic length 



(25) 



separating short length scales over which bending energy 
dominates from the large ones over which tension does. 

We now need to determine the force Ijl2|l on a parti- 
cle for the Hamiltonian described by Eqn. (|24|l . Using 
Eqns. HA3(1 and i|A4|l from the Appendix, we obtain 



f a = n(K ab - -Kg ab )K - ag ab e b - n{V a K) n (26) 



9S 



-K{Kl-Kf)-a 



(30) 



Note that the integrand has been decomposed with re- 
spect to a (right-handed) orthonormal basis adapted to 
the contour, {l,t,n}. 



C. Two-particle configurations 

We are interested in applying the general considera- 
tions of Sec. IIV Al to surface mediated interactions be- 
tween colloidal particles. In particular, we will consider a 
symmetrical configuration consisting of two identical par- 
ticles bound to an asymptotically flat surface, as sketched 
schematically in Fig. [21 

We label by {x, y, z} the Cartesian basis vectors of 
three-dimensional Euclidian space R 3 . Remote from the 
particles, the surface is parallel to the (x, y) plane. 

Let us agree that the constraining force fixes only the 
separation between the particles; their height, as well as 
their orientation with respect to the (x, y) plane are free 
to adjust and thus to equilibrate. This is also true of the 
contact line between surface and colloid when it is not 
pinned. Indeed, Kim et al. [l5| carefully argue that ver- 
tical forces and horizontal torques typically exceed hori- 
zontal forces and vertical torques by a significant amount. 
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FIG. 2: Two identical particles bound to an interface. As de- 
scribed in the text, the contour of integration can be deformed 
in order to take advantage of available symmetries. 




FIG. 3: Cross-section of a symmetric (solid line) and an an- 
tisymmetric (dotted line) two-particle geometry. 

Since the former can thus be assumed to very quickly 
equilibrate, they generally do not contribute to the mem- 
brane mediated interaction. 

There are two distinct manifestations of two-particle 
symmetry in this situation: either a mirror symmetry in 
the (y, z) plane (the symmetric case) or a twofold sym- 
metry axis, coinciding with the y axis (the antisymmetric 
case) . The former is relevant if the two particles adhere 
to the same side of the surface, the latter applies if they 
adhere on opposite sides (see Fig. [3J. In these two geo- 
metries the line joining corresponding points on the two 
particles lies along the z-direction. 

It is now possible to deform the contour of the line 
integral (|12fl to our advantage: as indicated in Fig. [3 the 
contour describing the force on the left hand particle may 
always be pulled open so that the surface is flat on three 
of its four branches (2, 3 and 4). The contributions from 
branch 2 will cancel that from 4; the only mathematically 
involved term stems from branch 1. The force on the 
particle is then given by 



dfi l a f a • 



(31) 



Let us now apply this general approach to a surface 
whose energetics can be described by the Hamiltonian 
density (|24"|) . 



D. The force between particles with symmetry 

1. Fluid membranes 

Both mirror and two-fold axial symmetry of branch 1 
imply that in Eqn. I|3(J|I the tangential term proportional 



to t vanishes. In the first case this follows from the fact 
that branch 1 becomes a line of curvature; hence, the 
curvature tensor is diagonal in (I, £)-coordinates and thus 
-ftTx|| vanishes. In the second case two-fold axial symme- 
try forces both K\\ as well as K± to be zero, since branch 
1 becomes a straight line and the profile is antisymmet- 
ric. In consequence, K = K± + K n = 0. We thereby 
obtain the first important simplification of the force from 
Eqn. I|30|) on that branch: 



n -n(Kl Kf) a 



I - k(V i K)n 



We now examine separately the two symmetric geome- 
tries (see discussion in Sec. II V C|) . 

a. Symmetric case. Tangent and normal vector on 
branch 1 lie in the (y, z) -plane, hence I = x. The deriva- 
tive of K in the direction of I along branch 1, W±K, is 
zero due to mirror symmetry. On branch 3 the surface is 
flat and thus the stress tensor is equal to f a 3 = —oe a . 
With this information we can calculate the total force 



Fx 



F S ymX on the particle: 



sym 



= aAL — Ik I ds (K 2 



(33) 



where AL > is the excess length of branch 1 compared 
to branch 3. If n = 0, we immediately have the im- 
portant general result that the force is always attractive 
irrespective of the detailed nature of the source. Unfor- 
tunately, the curvature contribution has no evident sign 
in general. However, for two parallel cylinders adhering 
to the same side of the interface the overall sign becomes 
obvious, as long as the particles are long enough such 
that end effects can be neglected: the contribution K 2 
then vanishes because branch 1 becomes a line. For the 
same reason AL = 0. This leads to the formula 



^sym,cyl/-^' 



2 J 



(34) 



where L is the length of one cylinder. Thus, the two 
cylinders repel each other. 

b. Antisymmetric case. Here branch 1 is a twofold 
symmetry axis and, as we have seen above, K\\ = K± = 
0. While the sign of Vx-fGi is not obvious, the deriva- 
tive V±K± is smaller than zero because K± changes sign 
from positive to negative. The profile on the midline 
is always tilted by the angle tp(s) in the direction indi- 
cated in Fig. |31 because any geometry with more than 
one nodal point in the height function between the par- 
ticles is expected to possess a higher energy. We fix the 
horizontal separation of the particles and allow other de- 
grees of freedom, such as height or tilt, to equilibrate (see 
Sec. II V C|l . The force on the particle is therefore parallel 



to x, F 



antisym 



F, 



antisym 



x : and given by 



F r 



antisym 



- r 



ct( cos <p(s) — l) 

nsmip(s) S7±(K ± + K {] )] , (35) 
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where we have used x ■ I — cos <p and x ■ n = — sin (p at the 
midpoint. Note that in this case the tension contribution 
is repulsive. As before, the sign of the curvature term is 
not obvious. 

If we restrict ourselves to the case of two parallel cylin- 
ders adhering to opposite sides of the interface, how- 
ever, then V_i_-ftT|| vanishes at the midpoint. Furthermore, 
|/ Q | is constant on each of the three free membrane seg- 
ments (due to Eqn. JBaJ). The stress tensor at branch 1, 
/ ~ f a ii must be horizontal to the x axis because ver- 
tical components equilibrate to zero as mentioned above. 
Let us look at the projection of the stress tensor onto I: 

h~f^-l={f X -x){x-l). (36) 

It follows that / • x = sign(/j/a; • I) \f \. We know that 
x ■ I = cosi/? > and fi = —a < 0. Hence, = — \f \x 
at the midpoint. This reduces Eqn. I|35|) to 

^a„tisym,cyl/£ = - O = y/ (T* + (kV ± K ± ) 2 - G > , 

which implies particle attraction. The length L is again 
the length of one cylinder. 

2. Membranes with tilt degree of freedom 

In Sec. IIIII we introduced a tangential vector field m a 
on the membrane, thereby modeling the degrees of free- 
dom associated with the tilt of the lipids. The minimal 
intrinsic Hamiltonian density Eqn. I|16|) already gives rise 
to a quite formidable additional metric stress, Eqn. Q17[l. 
Yet, for sufficiently symmetric situations the expression 
for the force simplifies quite dramatically, as we will now 
illustrate with another striking example. 

Let us consider two conical membrane inclusions which 
are inserted with the same orientation into a membrane 
at some fixed distance apart. Each inclusion will, due 
to its up-down-asymmetry, act as a local source of tilt. 
Provided the membrane is not in a spontaneously tilted 
phase, this tilt will decay with some characteristic decay 
length as described at the end of Sec. IIIII A typical situa- 
tion may then look like the one depicted in Fig. What 
can we say about the forces between the two inclusions 
mediated by the tilt field? 

Following the same reasoning as for the geometrical 
forces discussed above, and remembering that the tilt 
vanishes on branch 3 so that its contribution vanishes, 
we find 

F m = - J ds l a T^e b , (38) 

with T£ b given by Eqn. H17|) . To simplify this expression, 
we need to have a close look at the symmetry. For this 
it is very helpful to again expand vectors and tensors in 
a local orthonormal frame (I, t), just as we have done in 
the geometrical case above. Mirror symmetry then in- 
forms us that m" is an even function along the direction 




FIG. 4: Two conical inclusions act as sources of a local mem- 
brane tilt (inset). The tilt-field-lines are illustrated qualita- 
tively in this symmetric situation. 



perpendicular to branch 1, while m is an odd function 
and thus in particular zero everywhere on that branch. 
It thus follows that both Vj_77j" and vanish every- 

where on branch 1. Thus we have 

M = (Vim 1 ) + (Vym 11 ) , (39a) 
M ab M ab = (V^)^ (V||mll) 2 , (39b) 
e ab V a m b = Vxm 11 - Vp 1 = , (39c) 

where the "1" above the equation signs reminds us that 
this only holds on branch 1. We next need to look at 
the contractions of the individual terms in the metric 
material stress with l a e b . We find: 

L(V c m a )(V c m b )e b = {W±m x ) 2 l , (40a) 
L(V a m c )(V 6 m c )e b = (V ^fl , (40b) 
l a MM ab e h = M(V±m x )l . (40c) 

The two terms involving the derivatives m c V c can be 
rewritten by extracting a total derivative: 

l a m c (V c M)g ab e b = lm}V\\M 

= l\y\\{rJM)- (Vhto^M] .(41a) 

The total derivative will yield a boundary term once in- 
tegrated along branch 1, and since we assume that we are 
not in a spontaneously tilted phase, \m a \ will go to zero 
at infinity and thus the boundary term vanishes. With 
the same argument we find 

l a m c {V c M ab )e b = l[V\\(m ll (Vxm X )) 

-(Vnml^Vxm^)] . (41b) 

Again, the total derivative integrates to zero. Finally, 
the potential terms simplify to 

l a V(m 2 )g ab e b = V(m 2 )l , (42a) 
l a V'{m 2 )m a m b e b = . (42b) 
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Collecting all results, we arrive at the remarkably simple 
exact force expression F m — F m x, with 



(|A- 



els 



(V. 



(Viim 



J^ds V(r 



(43) 



There are two contributions to the force, one stemming 
from gradients of the tilt, the other from the tilt potential 
V. Remarkably, the tilt gradient contribution from each 
of the first two quadratic invariants has the same struc- 
tural form, thus the Lame coefficients A and /j, occur only 
as a combination in front of the integral. The modulus v 
has dropped out since the corresponding stress vanishes 
on the mid-curve (see Eqn. (|39c|) L 

The structural similarity of Eqn. (|43() to curvature 
mediated forces - Eqn. (|33|l - is very striking. Since 
| A + (J, > [3!|, the first integral states that perpendic- 
ular gradients of the perpendicular tilt lead to repulsion, 
while parallel gradients of the parallel tilt imply attrac- 
tions - the same "_L 2 — 1| 2 " motif as found in Eqn. (|33|l . 
Since in the untilted phase V(m 2 ) > 0, the second line 
shows that the integrated excess potential drives attrac- 
tion, just as the excess length (something like an inte- 
grated "surface tilt" ) drives attraction in Eqn. i.VMi . Un- 
fortunately, the overall sign of the force is not obvious. 
Looking at the field lines in Fig. the visual analogy 
with electrostatic interactions between like charged par- 
ticles would suggest a repulsion, but the above analysis 
advises caution (in Sec. IVlO we will see that this naive 
guess is at least borne out on the linearized level) . More- 
over, we should not forget that tilt does couple to geo- 
metry (namely, via the covariant derivative) and that the 
membrane by no means needs to be flat; hence, the con- 
tribution due to tension and bending given by Eqn. Q 
must be added, the sign of which is equally unclear. 



3. Further geometric Hamiltonians 

Within the framework of reparametrization invariant 
Hamiltonians providing a scalar energy density, a system- 
atic power series in terms of all available scalars and their 
covariant derivatives (each multiplying some phenomeno- 
logical "modulus" ) is a formal (and in fact standard) way 
of obtaining an energy expression of a physical system. In 
this respect the Hamiltonian (|24|l is no exception, being 
simply the quadratic expansion for an up-down symmet- 
ric surface (notice that a term proportional to K would 
break this symmetry, giving rise to a spontaneous curva- 
ture). We hasten to add that a second quadratic term, 
proportional to the Gaussian (or Ricci) curvature, exists 
as well, but this usually plays no role since it only results 
in a topological invariant (see also the Appendix). 

The fact that curvature (a "generalized strain") en- 
ters quadratically in the Hamiltonian i|24|) classifies this 
form of the bending energy as "linear curvature elas- 
ticity" (even though the resulting shape equations are 



highly nonlinear). However, for sufficiently strong bend- 
ing higher than quadratic terms will generally contribute 
to the energy densit y, g iving rise to genuinely nonlinear 
curvature elasticity [45(. Nevertheless, such effects pose 
no serious problem for the approach we have outlined so 
far. In fact, they are incorporated very naturally. We 
would like to illustrate this with two examples. 

a. Quartic curvature. Sticking with up-down symmet- 
ric surfaces, the next curvature order would be quartic, 
and this gives rise to three more scalars: K 4 , K 2 R and 
R 2 . Let us for simplicity only study the case of a quartic 
contribution of the form 



H 4 = \^iK 4 



(44) 



Using the general expression of the stress tensor for the 
scalar K n as calculated in the Appendix (see Eqn. ljA4|) ~) 
and going through the calculation from Scc. lIVDl we find 
for instance 



"'sym.cyl/-^ — 1_ ) 



(45) 



if two parallel cylinders adhere to the same side of the in- 
terface. This term increases the repulsion between cylin- 
ders found on the linear elastic level (see Eqn. I|34|| L pro- 
vided K4 > 0, i. e., provided the quartic term further 
stiffens the membrane. 

Assuming that 7i 4 perturbs the usual bending Hamil- 
we can use the two moduli to define a 



tonian ^kK 



characteristic length scale I4 := \/\k 4 \/k. The overall 
force up to quartic order can then be written as 



l±-{UK A _f 



(46) 



where the +- sign corresponds to stiffening. Notice that 
the correction term becomes only noticeable once the cur- 
vature radius of the membrane is no longer large com- 
pared to the length scale £4. It appears natural that £4 
is related to the membrane thickness, which for phos- 
pholipid bilayers is about 5 nm. Assuming a (quadratic 
order) bending stiffness of n ~ 20 k^T, we thus expect 
values for the modulus K4 on the order of 10 3 kBTmn 2 . 

b. Curvature gradients. In order length - it is possi- 
ble to also generate scalars which depend on derivatives 
of the surface curvature. One such term is 



H v = ^K V {V a K){V a K) 



(47) 



Using the expression for the stress tensor derived in the 
Appendix (see Eqn. (|A14|) ~) and again going through the 
calculation in Sec. IIVDI we find 



1 



1 



r v d^ 



(48) 



for the force between two symmetrically adhering cylin- 
ders. It depends on very subtle details of the membrane 
shape: the curvature is (roughly) a second derivative of 



10 



the membrane position, and this we need to square and 
differentiate two more times. Unfortunately, the sign of 
the interaction is not obvious here, as the second deriva- 
tive of Kj_ with respect to I may be either positive or neg- 
ative. Finally, we ca n also d efine a characteristic length 
scale here, £\r := \Z\ K v\/ K - The importance of a per- 
turbation 7iy of the usual bending Hamiltonian depends 
on whether or not the curvature changes significantly on 
length scales comparable to £y 



tensor identities. The right-hand side gives the numeri- 
cal values of the components of the covariant tensors gij 
and g lJ with respect to the coordinates x and y. 
The unit normal vector is equal to 



1 ( -Vh 
n = 

y/9 V 1 



(50) 



With the help of Eqn. the extrinsic curvature tensor 
is determined to be: 



V. DESCRIPTION OF THE SURFACE IN 
MONGE PARAMETRIZATION 

In the previous section analytical expressions for the 
force between two attached particles have been derived 
which link the force to the geometry of the surface at the 
midplane between them. It is worthwhile reemphasizing 
that they are exact, even in the nonlinear regime. In 
special cases, the sign of the interaction is also revealed. 

If one is interested in quantitative results, however, 
shape equations need to be solved - numerically or an- 
alytically. Either way, one needs to pick a surface 
parametrization. The choice followed in essentially all 
existing calculations in the literature is "Monge gauge" , 
and for analytical tractability its linearized version. The 
purpose of this and the following section is to translate 
the general covariant formalism developed so far into this 
more familiar language. To this end we first remind the 
reader what the basic geometric objects look like in this 
gauge. We arc then in a position to quantitatively study 
three different examples of interface mediated interac- 
tions in Sec. lVll 



A. Definition and properties 

Any surface free of "overhangs" can be described in 
terms of its height h(x, y) above some reference plane, 
which we take to be the {x, y) plane. Notice that x and y 
thus become the surface coordinates. The direction of the 
basis vectors {x, y, z} 6 K 3 is as described in Sec. IIVCI 

The tangent vectors on the surface are then given by 
e x = (1,0, h x ) T and e y — (0, l,h y ) T , where hi = dih 
S {x,y}). The metric is given by 



9ij 



h;hj 



(49) 



where Sij is the Kronecker symbol. Observe that is 
not diagonal; even though the coordinates {x, y} refer to 
an orthonormal coordinate system on the base plane, this 
property does not transfer to the surface they parame- 
terize. We also define the gradient operator in the base 
plane, V := (d x , d y ) T . The metric determinant can then 
be written as g = \gtj | = 1 + (Vft.) 2 , and the inverse met- 
ric is given by g lJ — 5ij — hihj/g. It is, perhaps, worth 
emphasizing that the latter, just as Eqn. I|49|) . are not 



V9 



(51) 



where — 



djdjh. 



Note that Eqn. I|51|) again is not a 
tensor equation; it provides the numerical values of the 
components of Kj in Monge gauge. 

Finally, it is also possible to write the trace K of the 
extrinsic curvature tensor in Monge parametrization: 



(52) 



B. Small gradient expansion 

In Sec. I VII we will be interested in surfaces that devi- 
ate only weakly from a flat plane. In this situation the 
gradient V/i is small, and it is sufficient to consider only 
the lowest nontrivial order of a small gradient expansion. 
K and dA can then be written as 



K 



-V 2 h 



0[{Vh) 2 ] , (53) 
dA = {l + i(V/i) 2 +C[(V/i) 4 ]} dxdy . (54) 



To evaluate the line integrals described in Sec. lIVDl we 
need expressions for K±_ and Kn as well as the derivatives 
V ±K± and Vx-^ii at branch 1 in Monge parametrization. 
In the small gradient expansion, the result is simply 



as well as 



K± 
Kn 



V±Kn 



-h xx {0,y) , 
-h yy (0,y) , 



-h 

-h 



x(o,y) 

yyx(0,y) 



(55a) 
(55b) 



(56a) 
(56b) 



We are now in a position to determine the forces between 
two particles in different situations. 



VI. EXAMPLES 

In this section we will illustrate the general framework 
of geometry-encoded forces by treating three important 
examples in Monge gauge: capillary, curvature mediated, 
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and tilt-induced interactions. For the first two, force- 
distance curves ha ve p reviously been derived on the lin- 
earized level 0, 0- EI ■ The route via the stress tensor 
reproduces these results with remarkable ease, thereby 
underscoring its efficiency and also confirming its va- 
lidity (at least on the linear level). To illustrate tilt- 
mediated interactions we restrict to a simplified situa- 
tion in which we neglect the coupling of membrane shape 
and tilt-order. Even if the geometry is "trivial" (a flat 
membrane), the material stress tensor is not, and forces 
remain. 

Both geometric examples are special cases of the 
Hamiltonian density (|24l) . When the gradients are small, 
the surface energy is given by the quadratic expression: 



2 



da; dy 



n(V 2 h) 2 + a{Vhf 



(57) 



If K — this describes a soap film; if k ^ it will describe 
a fluid membrane. 

The approach, traditionally followed in the literature, 
is to first determine the surface profile h(x, y) which min- 
imizes the energy Eqn. Q57J1 ■ For this one must solve the 
linear Euler-Lagrange equation 



V 2 (V 2 -r 2 )h(x,y) = Q , 



(58) 



where I is the length from Eqn. Ij25(l . In a next step, 
the energy corresponding to this shape is evaluated by 
reinserting the solution of Eqn. I|58(l into the functional 
Q57JI . This energy will depend on the relative positions 
of the bound objects. Appropriate derivatives of the en- 
ergy with respect to these positions will yield the forces 
between the particles. By contrast our approach - 
sidestepping the need to evaluate the energy — will be to 
determine the force directly from the surface profile using 
the line integral expressions for the force, Eqns. and 
(EH. 



A. Soap films 

For a soap film, n — and thus t = 0. The relevant 
Euler-Lagrange equation is therefore the Laplace equa- 
tion, V 2 h = 0. 

Consider first the symmetrical configuration consist- 
ing of two parallel cylindrical particles which adhere to 
one side of the soap film. Eqn. 133fl indicates that, if 
we neglect end effects, the force between the cylinders is 
proportional to the excess length on branch 1. The ex- 
cess length, however, is zero because the contact lines are 
straight. Therefore, the force is also zero. Likewise, in 
the antisymmetric configuration with adhesion on oppo- 
site sides, the soap film between the cylinders will be flat 
if the vertical particle displacements are allowed to equi- 
librate. Therefore, (p(s) = (see Fig. |3J) and Eqn. J35J) 
will yield a zero force exactly as in the symmetric case. 
In an analogous way one obtains the same result for the 
case of two spheres. 




FIG. 5: Definition of the coordinates for a single quadrupole 
(viewed from above). 



The situation is less simple if the film is pinned to the 
particle surface. Let us consider two spherical particles 
of radius a with a contact line that departs only weakly 
from a circle. 

Stamou et al. Q have studied this case by using a su- 
perposition ansatz in the spirit of Nicolson 46] : first, 
the height function of one isolated particle is determined 
with the correct boundary conditions. Then, the com- 
plete height function is assumed to be the sum of the two 
single-particle fields of each of the two colloids. Strictly 
speaking, this approach destroys the boundary conditions 
at the particles' contact lines; it does, however, give the 
correct leading order result for large separation |47|. 

Using polar coordinates p and 0, the solution of the 
shape equation outside a single spherical particle can be 
written as @ 

ft S phere(p, <t>) = ^° llX (~) 
oo 

+ } j A m cos [m(<j> - (p m ,o)] ( - ) ,(59) 

m— 1 

with multipole coefficients A m and phase angles 4> m ,o- 
The former can be determined as follows: The monopole 
Ao vanishes because there is no external force such as 
gravity pulling on the particle. The dipole coefficient A\ 
characterizes the tilt of the contact line relative to the z 
axis; it also vanishes if there is no external torque acting 
on the sphere. All higher multipole coefficients can be 
read off from the Fourier expansion of the contact line 
at p = a. It is intuitively obvious and indeed confirmed 
by a more careful calculation 0, @ that the quadrupole 
dominates the energy at lowest order. 

One can therefore restrict the calculation to the single- 
particle height function 

/ a \ ^ 

fr s phere(P,<£) = Qcos[2(0 - (f>o)}{-) , (60) 

where </>o :— 02. o is the angle that represents the rotation 
of the particle about z (see Fig. [SJ. 

If the complete height function is a superposition, as 
described above, the force on the left particle in lowest 
order has been found to be 
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FIG. 6: Two quadrupoles on a soap film (viewed from above). 



for the symmetric (<fo,A = — <Po.b) and the antisymmetric 
(0o, A = 0, 0o, b = tt/2) configurations (see Fig.^J. 

Let us now examine the same two configurations using 
the line integral representation for the force. 

a. Symmetric case. The force Eqn. I|33[l is proportional 
to the excess length of branch 1 with respect to branch 
3. To quadratic order in gradients, this difference can be 
written as 



AL = lim 



lim 



L/2 

-L/2 
L/2 

-L/2 



Ay 



Ay 



l + Ag(0, tf )-l]} 



-h 2 y {^y)+0[{Vhf]]} 



(62) 



The height function along the symmetry line between the 
particles can be expressed in Cartesian coordinates as: 

r,2 



h(0,y) = 2Qcos 



2( arctan 



2// 



</>0,A 



y 



4 



(63) 



Substituting into the second line of Eqn. H62J) gives 

L 



AL = 



lim 

L — >oo 

48ttQ 5 



96Q Z 



d 5 



— arctan - , 
d 5 d 



0(L~ 



(64) 



implying via F t 



sym,soap 



= <jAL the same force as obtained 
from energy minimization, Eqn. I|61[) . However, it would 
be fair to say that we have gained additional information 
concerning the nature of this force, missing before. The 
force is directly proportional to the length added to the 
mid-curve as it is stretched. A geometrical interpretation 
has been provided for the force. Recall also that this is a 
non-perturbative result: it does not depend on the small 
gradient approximation. 

b. Antisymmetric case. In this case the horizontal 
force is given by Eqn. (|35|l with k = 0: 



antisym.soap 



cr lim 

L — >oc 



a lim 

L — >oc 



a lim 

L — >oo 



r f L/2 






n ■ z — 1 


L J -L/2 





L/2 

-L/2 
L/2 

-L/2 



Ay 



dy 



} 



-lhl(0,y)+O[(Vh) i }]} .(65) 



The height function between the particles is given by 
h(x, y) = Qa 2 x 



cos [2 (arctan 

y 2 + (l + x) 2 



)] cos [2(arctan -rM] 

2 

2/ 2 + (f -x) 2 



Mo, i/) 



32Qa 2 d(d 2 - 12y 2 ) 



(d 2 + Ay 



2\3 



(66) 



(67) 



Inserting this into Eqn. I|65l) yields a force which again 
agrees with the one obtained in Eqn. I|61l) . 

As an example, let us look at colloids with a radius of 
1 /xm trapped at the air- water interface (a ~ 70mN/m), 
which have a pinning quadrupole of 1% of their radius 
(Q ~ 10 nm). At a separation of 3/Ltm they feel an (at- 
tractive or repulsive) force of 1 pN, and at a separation 
of about 16 /xm their interaction energy is comparable 
to the thermal energy. These forces are not particularly 
strong, but they act over an exceptionally long range. 



B. Fluid Membranes 

To describe a fluid membrane, it is necessary to in- 
clude the bending energy in Eqn. l(57|) . Let us focus on 
the problem of two parallel adhering cylinders which are 
sufficiently long so that end effects can be neglected (the 
fluid membrane analogue of the problem examined for 
soap films). In this case the height function of the sur- 
face depends only on one variable, x. Recall that for the 
corresponding soap film case no interaction occurred (in 
the absence of pinning), see Sec. I VI Al 

a. Symmetric case. Using the energy route, Weikl 
[ll| shows that, at lowest order in the small gradient 
expansion, the energy per unit length of the cylinder is 

m 



(K + 2i? 2 C/) 2 (tanh^ - 1 



(68) 



Here R is the cylinder radius, U is the adhesion energy 
per unit area, I is the characteristic length defined in 
Eqn. 1)25(1 . and d is the distance between the two centers 
of the cylinders. To obtain the force per unit length L of 
the left cylinder, we differentiate Eqn. I(68|l with respect 
to d Ei: 



1 / n + 2R 2 U 



-^ l/L ~ 2 K V2 Kj Rcosh^ 



(69) 



The cylinders always repel. 

We would now like to determine the force using the line 
integral of the corresponding stress tensor. Rewriting the 
relevant Eqn. I|34|) in small gradient expansion yields (see 
Eqn. JSSaJ): 



i,cyl/i 



1 



^L(o) 



(70) 
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We use the expression for h given in Ref. [Tl|: 

, N (K + 2i? 2 [f)coshf , s 

h(x) = — L + const . 71 

2cri?cosh^ 

Its second derivative with respect to x at x = is 

K + 2R 2 U 



h xx (0) 



2«;i?cosh ^ 



(72) 



Inserting this result into Eqn. I|7U|I reproduces the force 
given by Eqn. I|t)9[) . 

b. Antisymmetric case. For two cylinders on opposite 
sides of the membrane the energy is given by 0, \^ 



E; 



antisym 



,cyl(d) 



(« + 2i? 2 C/) 2 (coth^ - 1 
IJ^RR 2 



which gives a force on the left cylinder |4£ 



antisym, cyl 



/ L = -K 



1 / k + 2R 2 U x 2 



2 V2Ki?sinh 



2f 



(73) 



(74) 



The small gradient expansion of Eqn. 13 7|) is 

F.. 



„.>,:, is ym,c y i/L = ^K(eX7 ± Kj_) 2 = -n[£h xxx (0)f 



If we now again take the height function from Ref. [1 
we arrive at 



h(x) 



{k + 2R 2 U) sinh § 
2cri?smh^ 



which yields 



£^xxx (0) 



2R 2 U 



2«i?sinh^ 



(76) 



(77) 



Inserting this into Eqn. (|75(1 reproduces the result ()74() . 

How big are these forces? As an example, let us look 
at an actin filament (R ~ 4nm) adsorbed onto a mem- 
brane with a typical bending stiffness k ~ 20 fceT, where 
k^T is the thermal energy. Noting that \J2XJ j k will be 
the contact curvature at the point where the membrane 
detaches from the adsorbing filament [5(j and that this 
should not be too much smaller than the bilayer thick- 
ness in order for a Helfrich treatment to be permissible, 
we take 2UR 2 /k ~ 1 as an upper limit. We then find that 
two adsorbed actin filaments at a distance d ~ I (where 
approximately sinh w cosh ps 1) feel a force of about 
2 — 3pN/nm. Alternatively, we can calculate at what 
distance the interaction energy per persistence length of 
the filament (£ p w 15 /im) is of order k^T. Using a typ- 
ical value for cell membranes of I ~ 50 nm, we obtain a 
separation of about 0.7 /im. This is huge, and should re- 
mind us of the fact that on this scale a lot of membrane 
fluctuations will occur which we have neglected. Still, it 
shows rather vividly that membrane mediated forces can 
be very significant. 



C. Lipid tilt 



The discussion in Sec. IIIII shows that lipid tilt order, 
described by the surface vector field m a , influences the 
shape of the membrane, even if the Hamiltonian den- 
sity does not contain an explicit coupling of m a to the 
extrinsic curvature. The coupled system of differential 
equations (|18|l poses a formidable task, clearly exceeding 
the already substantial one for the undecorated shape 
equation alone. 

Our priority is to illustrate the workings of the gen- 
eral formalism, therefore we will limit the discussion to 
a simple case where the analytical treatment is rather 
transparent: we will assume that the membrane itself re- 
mains flat, such that the energy density stems exclusively 
from lipid tilt (as described by H m from Eqn. (ftT))l ). This 
is not a self-consistent approximation, but should give a 
good description in the limit in which the tilt moduli A 
and fj, are significantly "softer" than the bending modu- 
lus. In this case the inclusions we have talked about in 
Sec. II V lT2l will predominantly excite tilt and not bend. 
More sophisticated (analytical and numerical) studies of 
lipid tilt and mediated interactions exist, which provide 
better quantitative answers |24| . 

For flat membranes, the Euler-Lagrange equation itT^l) 
reduces to 



(A + Ai)VV • m - 



fiV 2 m 



2V'm = 



(78) 



where m is the 2d tilt vector in the membrane plane. Fo- 
cusing first on one inclusion, the situation acquires cylin- 
drical symmetry. Writing m(r) = m(r)e r and restricting 
to the untilted membrane phase, for which the tilt po- 
tential is sufficiently well represented by V{m 2 ) = ^tm 2 
with t > 0, Eqn. I|78() reduces to a simple Bessel equation 



2 /. 
x m 



(x 2 



Dm = 



(79) 



where x = r/£ m , £ m is the length defined in Eqn. J22J>, 
and the prime denotes a derivative with respect to x. The 
solution is 



m(r) = mo 



ifi(r /C) 



(80) 



where tq is the radius of the inclusion, mo the value of 
the tilt at this point, and K v a modified Bessel function 
of the second kind 51]. As anticipated, the tilt decays 
essentially exponentially with a decay length of £ m . 

Obtaining the exact tilt field for two inclusions is very 
difficult, since satisfying the boundary conditions is trou- 
blesome. However, if we again use the Nicolson approx- 
imation |46j | and assume that the total tilt distribution 
is given by the superposition of two solutions of the kind 
(|80|l . things become manageable. The tilt-mediated force 
between two symmetric inclusions is then obtained by in- 
serting the appropriate values and derivatives of the tilt 
field m(x,y) on the mid-line into Eqn. 14311 . After some 
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straightforward calculations we get the force 52] 

f°° 1 
F m = At£ m m* 2 / d£ — - x 

Jd- we - d * 2 

\e - 2d* 2 )K (OK 2 (o + (e - d* 2 )K 2 (o 

= -1-KtUm^K x {dlt m ) . (81) 

where = mo/Ki(ro/£ m )i d* — d/2£ m , and d is the 
separation between the inclusions. As we see, the force is 
repulsive and decays essentially exponentially with dis- 
tance over a decay length of £ m . Integrating it, we get 
the repulsive interaction potential 

U m (d) = 2nU 2 a m* 2 K Q (d/£ m ) . (82) 

Let us try to make a very rough estimate of how big 
such a force might be. For this we need to obtain some 
plausible values for the numbers entering into Eqn. (|81|l . 
For t we may use the equipartition theorem and argue 
that ^t(m 2 )a = ^kftT, where a is the area per lipid and 
k^T the thermal energy. Assuming that the root-mean- 
square fluctuations of m are 10° and using the typical 
value a ~ 0.75 nm 2 , we get t ~ 40/cB7 1 /nm 2 . Assum- 
ing further a rather conservative tilt decay length of the 
order of the bilayer thickness, i. e. £ m ~ 5 nm, that the 
inclusion has a radius of ro — 3 nm and imposes there 
a local tilt of m ~ 0.2, we find that two inclusions at a 
distance of 10 nm feel a significant force of about 17pN. 
And at a distance of d ~ 22 nm their mutual potential 
energy is 1 k^T compared to the separated state. Notice 
that this is much larger than the Debye length in phys- 
iological solution, which is typically only 1 nm. Hence, 
tilt-mediated forces can compete with more conventional 
forces, such as (screened) electrostatic interactions. It 
should be kept in mind, however, that if we permit the 
membrane to bend, some of the tilt strain can be relaxed, 
thereby lowering the energy. 



VII. CONCLUSIONS 

We have shown how the stress tensor can be used to 
relate the forces between particles bound to an inter- 
face directly to the interface geometry. In this approach, 
the force on a particle is given by a line integral of the 
stress along any closed contour surrounding the particle. 
The stress depends only on the local geometry; thus the 
force is completely encoded in the surface geometry in 
the neighborhood of the curve. 

The relationship between the force and the geome- 
try provided by the line integral is exact. In the lin- 
ear regime, as we have shown for selected examples in 
the previous section, the force determined by evaluating 
this line integral reproduces the result obtained by the 
more familiar energy based approach. Unlike the latter, 
however, our approach permits us to consider large de- 
formations. The expression for the line-integral is fully 



covariant, involving geometrical tensors; one is not lim- 
ited to any one particular paramctrization of the surface 
such as the Monge gauge. Indeed, as we have seen the 
geometrical origins of the force can get lost in this gauge. 

As we have emphasized previously, this approach is not 
a substitute for solving the nonlinear field equations. To 
extract numbers, we do need to solve these equations. 
But even before this is done, the line integral expression 
can provide valuable qualitative information concerning 
the nature of the interactions between particles. This is 
because the geometry along the contour is often insensi- 
tive to the precise conditions binding the particle to the 
interface. This contrasts sharply with the energy based 
approach; there, one needs to know the entire distribu- 
tion of energy on the interface before one can say any- 
thing about the nature of the interaction. As we have 
seen in the context of a symmetrical two-particle config- 
uration, it is sometimes relatively easy to identify quali- 
tative properties of the geometry] it is virtually impossi- 
ble to make corresponding statements about the energy 
outside the linear regime. 

The stress tensor approach also has the virtue of com- 
bining seamlessly with any approach we choose, be it 
analytical or numerical, to determine the surface shape. 
Thus, for instance, one can find surfaces that minimize 
a prescribed surface energy functional using the program 
"Surface Evolver" |5^| • The evaluation of the force via a 
line integral involving the geometry along the contour is 
straightforward; in contrast, the evaluation of the energy 
involves a surface integral, and the forces then follow by 
a subsequent numerical differentiation. In other words, 
the route via the energy requires one more integration 
but also one more differentiation. This appears neither 
economical nor numerically robust. 

We have illustrated how internal degrees of freedom on 
the membrane can be incorporated within this approach 
using a vector order parameter describing lipid tilt as an 
example. It is indeed remarkable just how readily non- 
geometrical degrees of freedom can be accommodated 
within this geometrical framework. Here again, new ex- 
act non-linear expressions for the force between particles 
mediated by the tilt are obtained which are beyond the 
scope of the traditional approach to the problem. Vari- 
ous patterns emerge which could not have been guessed 
from inspection of the Hamiltonian, in particular the ex- 
istence of a "_L 2 — 1| 2 " motif common to the geometrical 
and tilt mediated forces between symmetrical particles. 

We have considered the force between a pair of par- 
ticles. However, the interaction between more than two 
particles is generally not expressible as a sum over pair- 
wise interactions; superposition does not hold if the the- 
ory is nonlinear (see Ref. ^5| for a striking illustration). 
This, however, poses no difficulty for the stress tensor 
approach, because the underlying relation between sur- 
face geometry and force is independent of whether or not 
a pair-decomposition is possible (see Fig. 0). For cer- 
tain symmetric situations a clever choice of the contour 
of integration may again yield expressions for the force 
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APPENDIX A: 




FIG. 7: Three-body interactions. The force on one particle 
can be obtained by integrating the surface stress tensor along 
a line of integration enclosing that particle, (cf. Eqn. 112H '). 



analogous to those obtained in Sec. 11V 131 

Multi-body effects become particularly relevant when 
one considers 2D bulk phases as, for instance, in a system 
consisting of a large number of mutually repulsive parti- 
cles adhering to one side of an interface. In this case it 
is possible to identify expressions for state variables such 
as the lateral pressure by exploiting the approach which 
has been introduced here. 

The interfaces we have considered are asymptotically 
flat and thus support no pressure difference. At first 
glance it may appear that our approach fails if there is 
pressure because the stress tensor is no longer divergence 
free on the free surface (one has V a / a = Pn) which 
would obstruct the deformation of the contour described 
in Sec. IIV CI As we will show in a forthcoming publi- 
cation, however, it is possible to adapt our approach to 
accommodate such a situation. 

The interactions we have examined correspond to par- 
ticles whose orientations are in equilibrium. The ad- 
ditional application of an external torque {e.g. on two 
dipoles via a magnetic field) will introduce a bending mo- 
ment. It is, however, possible to treat such a situation 
with the tools provided in Ref. : Just as translational 
invariance gives us the stress tensor, rotational invariancc 
gives us a torque tensor. Its contour integrals provide us 
with the torque acting on the patch one encircles. 

Finally, genuine capillary forces involve gravity. The 
associated energy, however, depends not only on the geo- 
metry of the surface but also on that of the bulk (it is a 
volume force). The results thus differ qualitatively from 
those presented here. This will be the subject of future 
work. 
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In Sec. Ill Dl we derived the general expression 

f a = (T ab - Ti ac K b )e b - V b Ti ab n (Al) 

for the surface stress tensor. In this appendix we will 
specialize <|A1|) to a few important standard cases. 

a. Area. The simplest case is the area, Ti = 1, which 
is (up to a constant prefactor) the Hamiltonian density 
of a soap film. We evaluate Ti ab and T ab appearing in 
Eqn. (53) using Eqn. 0: H ab = 5Ti/5K ab = and 



-g ; for T ab we use the identity 



Thus we get 



IT- = \^ 

og a b 2 v 



f a = -9 ab e b 



(A2) 



(A3) 



Note that the functional derivatives S in this case are 
equal to partial derivatives d because Ti does not depend 
on higher derivatives of g ab or K ab . 

b. Powers of K . For the Hamiltonian density Ti = 
K n = ( g abK ab ) n one derives 54]: H ab = nK n 
T ab = 2nK n - 1 K ab - K n g ab which gives 



1 g ab and 



f a = {nK n - 1 K ab - K n g ab )e b - n(V a K n - l )n . (A4) 

The case n = 2 is needed in Eqn. (|24Jl . 

c. Einstein- Hilbert action. Canham [42j originally 
used the quadratic Hamiltonian Ti = K ab K ab . For 
this one we easily see that Ti ab = 2K ab and T ab = 
-Tig ab + 4K^K bc [54|. Using the contractions of both 
Gauss-Codazzi-Mainardi equations © as well as the 
fact that the Euler-Lagrange derivative £(Ti) is linear 
in the Hamiltonian, we get with the help of Eqn. © the 
Euler-Lagrange derivative of the Einstein-Hilbert action, 
Ti = R: 



£(R) = £(K 2 ) - £(K ab K ab ) = -2K ab G ab 



(A5) 



Here, G ab = R ab — \Rg a b is the Einstein tensor, which 
vanishes identically in two dimensions. Thus, surface 
variations of K 2 and K ab K ab differ only by boundary 
terms (in accord with the Gauss-Bonnet theorem 25]). 
In higher dimensions, however, £(R) oc G ab is a non- 
trivial result, and the above seemingly inelegant (since 
extrinsic) derivation is after all remarkably economical. 

d. Curvature gradient. The next example we con- 
sider is the Hamiltonian density Ti = |(V C K)(V C K) = 
i(VA') 2 . Now we need to keep in mind that Ti ab and 
T ab are functional derivatives 



Ti 



ab 



5Ti 



&H 



SKab dK ab 



dTi 



dV c K al 



(A6) 



because Ti depends on derivatives of K ab . We obtain 

Ti ab = -X7 c (g ab X7 c K) = -g ab AK . (A7) 
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The determination of T ab is a little more difficult; to 
avoid errors, let us proceed cautiously and consider the 
variation of the Hamiltonian H = h J dA (X7K) 2 with 
respect to the metric g ab and identify T ab at the end of 
the calculation. The variation yields: 

5 g H = \J d 2 £S g [^g(VK) 2 } (A8) 

= \jdA{^L{VKY+5 g [{VKf]}. 

The evaluation of the first term involves the reuse of 
Eqn. 1|A2(1 : we expand the second term (5f| 



6 g [(WK) 2 } = 5 g [g ab ](V a K)(V b K) 

+2g ab {V a K)5 g [V b {K cd g cd )] 

= -{V a K){V b K)8g ab 

+2g ab (V a K)\7 b (K cd 5g cd ) 

= -(V a K)(V b K)6g ab 

-2g ab (V a K)V b (K cd Sg cd ) , (A9) 

where the identity 5 g g ab = —g ac g bd Sg cd is exploited 
twice. We obtain 



S g H = i / ,1.1 



^g ab (VK) 2 ~(V a K)(V b K) ,*,,„,, 



dA g ab (X7 a K)X7 b (K cd 5g cd ) . (A10) 
The last term can be rewritten as 

J dA g ab {V a K)V b {K cd 5g cd ) = - J dA K ab AKSg ab 



dA V f 



g ab {V a K)K cd 5g cd 



(All) 



the second term on the right hand side is a total diver- 
gence and can be cast as a boundary term. Therefore, it 
does not contribute to T ab . Collecting results, we find 

S g H IP -1 J dA T ab 5g ab + boundary terms , (A12) 

with 

T ab = (ya K j (yb K ^ _ i^VJf) 2 - 2K ab AK . (A13) 

Thus, for H = i(Vif) 2 we get for f a given by Eqn. l|AT)l 
the remarkably compact expression 



(V a K)(y b K) - ^g ab {S7K) 2 ~ K ab AK e b 
- V a AK n . (A14) 



e. Vector field. As a final example let us consider 
Hamiltonians of the kind introduced in Sec. IIIII which 
have internal vector degrees of freedom. With the sym- 
metric tilt-strain tensor M ab = i(V a m b + V b m a ) we can 



for instance look at the quadratic Hamiltonian density 
H = i(V a m a ) 2 = \M 2 , where m a is the (contravari- 
ant) surface vector field and M = g ab M ab . This term 
is purely intrinsic, hence Ti ab = 0. The difficult part is 
the covariant differentiation, which acts on a vector field 
and is thus dressed with an additional Christoffel sym- 
bol. Since the latter depends on the metric and its first 
partial derivative |26| , it will contribute to the variation: 



dA 



5 g K a + n b m b ) 2 }(Al5) 



The first term is once more simplified via Eqn. ( | A2 
while the second term calls for the Palatini identity Q 



°g l ab 



1 



g cd [V b Sg da + V a 6g bd - V d 5g ab ] 



(A16) 



Since the 5g ab are the components of a tensor (they must 
describe a proper variation of the metric tensor), the vari- 
ation S g T ab is also a tensor, even though the Christoffel 
symbol itself is not. Using Eqn. (|A16fl . the second term 
m Eqn. (|AT5)l can thus be rewritten as 



b\2 



6 g (m° a + T a ab m») 



2Mm b S g T a ab 



Mm b g ad {V b 6g da + V a Sg bd - V d 5g ab ) 



Mm d g ab (V d Sg ab ) . 



(A17) 



The derivative of Sg ab is removed by a final partial inte- 
gration. Collecting everything, we thus find (up to irrel- 
evant boundary terms) 



S g H 



dA 



{Ul 2 -V d [m d M]}g ab Sg ab 
{ - \m 2 - m d V d M }g ab 6g ab 



Thus, the metric stress tensor is 



rjiab _|_ 



M 2 + 2m c V c M 



ab 



(A18) 



(A19) 



Notice that it is directly proportional to the metric; its 
effect in the stress tensor will thus be to renormalize the 
surface tension. 

The second quadratic invariant, H = M ab Al ab , does 
not provide any additional difficulties compared to \M 2 , 
even though the calculation is a bit longer. One finds: 

T ab = -M cd M cd g ab + 2MM ab + 2m c V c M afc 

- (V c m Q )(V c m h ) + (V a m c )(V h m c ) . (A20) 

Finally, the third quadratic invariant TL = jF ab F ab 
(with F ab — V a m b — V h m a ) can be treated rather easily 
by noting that F ab = d a m b — d b m a is independent of the 
connection. A short calculation then shows that 



T 



ab 



pacpb 



1 



g ab F cd F 



cd 



(A21) 
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This is nothing but the energy-momentum tensor from 
electrodynamics 36] . In two dimensions it can be further 
simplified, since any antisymmetric tensor is then propor- 
tional to the epsilon-tensor: F ab — he EdF . Inserting 
this into Eqn. I|A21|I and using the identity e ac e b c = g ab 
|37j . we find 

T ab = ^ gab ^ £cdX7Cmd y ^ (A22) 



showing that the stress is isotropic, just as in the case of 
the Hamiltonian Ti — \M 2 . It will thus only renormalizc 
the surface tension and, in particular, not single out any 
specific new directions on the membrane. 
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